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^ ■ ABSTRACT 
O \ 

^ \ We study non-LTE inversions of the Ca II infrared triplet lines as a tool for 

, inferring physical properties of the quiet Sun. The inversion code is successful 

\ in recovering the temperature, velocity and longitudinal magnetic flux density 

^ \ in the photosphere and chromosphere, but the height range where the inversions 

O ! are sensitive is limited, especially in the chromosphere. We present results of 

^ \ inverting spectropolarimetric observations of the lines in a quiet Sun region. We 

^ \ \ find three distinct ranges in chromospheric temperature: low temperatures in the 

^ \ internetwork, high temperatures in the enhanced magnetic network and interme- 

\ diate temperatures associated with low magnetic flux regions in the network. 

^ \ The differences between these regions become more pronounced with height as 

\ the plasma-/? decreases. These inversions support the picture of the chromo- 

^ \ sphere, especially close to the magnetic network, being highly inhomogeneous 

\ both in the vertical and horizontal directions. 
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Introduction 



A transition from a high plasma-/^ (/? = Pgas/Pmag) to a low-/? regime takes place in 
the chromosphere. With this transition the nature of the dynamics changes: in the high-/? 
regime the plasma dominates the magnetic fields whereas in the low-/3 regime the opposite 
is true. Waves undergo mode mixing, reflection and refraction in the layer where plasma- 
(3^1. Because of the tight connection between the magnetic fields and dynamics it is 
crucial to understand how the two are coupled in t he quiet Sun chromosphere. Several 
recent observational works have focused on this (e.g., 



Judge et al 



2001 



Mcintosh &: Judge 



200ll : iMcIntosh et all l2003l : Ijefferies et all 12009 and IVecchio et all 120071 ). All find a clear 
correlation between the magnetic field topology and waves in the chromosphere. None of the 
above mentioned works, however, directly addresses the magnetic field in the chromosphere. 
Instead, proxies such as potential field extrapolation or the existence of fibril-like structures 
in the emission of chromospheric lines, are used. 

In order to obtain physical parameters (instead of proxies) from spectropolarimetric 
observations an inversion problem is unavoidable. Since the radiative transfer equation 
for polarized light is a set of coupled differential equations, and in the chromosphere local 
thermodynamical equilibrium (LTE) is no longer valid, solving the inverse problem is chal- 
lenging. Furthermore, one is never certain that any solution obtained is in fact the correct 
solution, instead we have to rely on some a priori knowledge of the solar atmosphere. Inver- 
sions are commonly used to to produce quantitative information of the photosphere. They 
are not yet as widely used for chromospheric observations. Significant pr ogress in chromo- 



spheric inversions has been made for the He I 10830 A multiplet ( e .g., Lagg et al 



Centeno et al. 


2005; 


Trujillo Bueno et al. 


2005; 


Sasso & Solanki 


2006: 


and Ca II infrared 1 


IR) triplet lines (e.g.. 


Socas-Navarro et al. 


2000a 



Socas-Navarro et al.ll2006l ). 



20041 : 

bood) 

Socas-Navarrq i2005i : 



The magnetic structure of the quiet Sun chromosphere is still a matter of debate, but 
there is no doubt that it plays an important role. The leakage of low- frequency (i.e., below 
acoustic cut-off frequency) waves into the chromosphere and the formation o f spicules have 
been attributed to the existence of inclined magnetic fields in the network (e.g. , iDe Pontieu et al, 
20041 : Ijefferies et al.ll2006l : iHansteen et al.ll2006l : iDe Pontieu et al.ll2007). It is also suggeste d 
that these waves may be related to the heating of the chromosphere ( Ijefferies et al.ll2006l ). 
Existence of canopy-like structure s such as proposed by lGabriell fll976l ) is still an open debate. 
For example IVecchio et al.l (120071 ) observe oscillatory patterns in t he chromosphere t h at are 
hig hly consistent with the c anopy- scenario, while models by e.g., ISchrijver fc Title! (120031 ) 
and IJendersie fc Peterl (120061 ) show that the existence of a canopy-like structure is disrupted 
by internetwork magnetic fields. The connection between magnetic fields, dynamics and 
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thermal properties of the quiet Sun chromosphere appears complex. In this paper we extend 
the use of non-LTE inversions to the quiet Sun chromosphere. With non-LTE inversions we 
do not have to rely on proxies, such as line intensity, to infer physical quantities. 

This paper is a continuation to Pietarila et al. (2006, hereafter P06) and Pietarila 
et al. (2007, hereafter P07) in which the use of Ca II IR triplet lines as a diagnostic for 
chromospheric magnetism are studied. Here we focus on non-LTE inversions of the Ca II IR 
triplet lines, and what information can be extracted from them. In section [2] the non-LTE 
inversion code used in this work is described in some detail. Also the sensitivity of the 
inversions is discussed. We demonstrate that the Ca lines are sensitive only in a narrow 
height range in the chromosphere, and that our inversion code cannot recover discontinuities 
in magnetic field inclination or velocity (at least with the observables used). We also show 
that discontinuities can cause complex patterns in Ca line Stokes V asymmetries. Inversions 
of quiet Sun observations are discussed in section [31 The network is found to be more dynamic 
than the internetwork. Three different regimes are found in chromospheric temperature: low 
temperatures corresponding to the "magnetic field free " internetwork, high temperatures 
in the enhanced magnetic network and an intermediate temperature regime associated with 
low magnetic flux regions. The main results are discussed (in the context of previous work) 
in section HI and final remarks are made in section [51 



2. Inversion code 



The n on-LTE inversion code employ ed here has been described in detail in lSocas-Navarro et al 



( 1l998l ) and lSocas-Navarro et al.l (l2000bl ). For completeness, we provide here a brief summary 
of its main features. The user needs to supply an initial model atmosphere, which consists 
of the run with standard optical depth (r5oo) of temperature, microturbulence, line-of-sight 
velocity and magnetic field vector (strength, inclination and azimuth). In addition to these 
vector quantities, a model atmosphere includes the following single-valued parameters: elec- 
tron pressure at the top of the atmosphere, macroturbulence and filling factor. The macro- 
turbulence is used also to describe the spectral instrument profile. The filling factor is used 
to consider that the magnetic structure is smaller than the resolution element (or that there 
is some stray light in the spectrograph). 

With the top boundary condition for the electron pressure and the run of tempera- 
ture, the hydrostatic equilibrium equation is solved to obtain the run of electron pressure, 
gas pressure and density and the geometrical scale height. The code then uses a precondi- 
tioning strategy (iRybicki Sz Hummerlll99ll : ISocas-Navarro fc Trujillo Buend 119971 ) to solve 
the statistical equilibrium equations in that atmosphere and synthesizes the Stokes profiles 
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emerging from the initial model as well as the response functions, i.e. the derivatives of the 
Stokes vector at a given wavelength with respect to the model parameters. 

The derivatives enter a Levenberg-Marquardt scheme that determines the first-order 
correction needed to compensate for the misfit between the observed and the synthetic pro- 
files. The correction for a given depth-dependent parameter (e.g., temperature) is obtained 
as a smooth function of t^qq, depending on the number of nodes selected by the user for that 
parameter (for two nodes the code retrieves a correction that is linear with t^qo, for three 
nodes a parabola, and for four or more nodes the correction is a cubic spline interpolation). 
The entire algorithm is iterated until it converges to a solution that minimizes the difference 
between the observed and synthetic profiles. 

2.1. Sensitivity of the inversion code 

The inversion code retrieves model atmospheres with the atmospheric parameters as 
a function of optical depth, ranging from log^r) = [—6, 0.8]. The sensitivity of the spec- 
tral lines, however, is limited to only a range of atmospheric depths. To identify what the 
depths are we made the following numerical experiment: a set of simultaneous Stokes I 
and V profiles of the 8498 A and 8542 A Ca II lines and the two photospheric Fe I lines, 
8497 A and 8538 A (located in the wings of the Ca lines), were inverted using 100 different 
initializations. The initial models were chosen randomly from existing inversions of observa- 
tions. To make the test comparable to the observations discussed later, we do not include 
Stokes Q 01 U profiles. Additionally, we force the magnetic field to be vertical, i.e., only the 
longitudinal magnetic field is considered. The lines are, as expected, sensitive to only part 
of the atmosphere (fig. [T]). The spread in the output models is smaller in two regions: in 
the photosphere, where the Fe lines and the Ca line wings are formed, and in the chromo- 
sphere, where the Ca line cores are formed. For temperature the heights are log{T) and 
log{T) ^ —5 , for magnetic flux log{T) ~ —1 and log{T) ~ —5 , and for velocity log{T) ^ —1 
and log{T) ^ —5.5. Since neither the Ca nor the Fe lines are formed in the region between 
logij) ~ [—4, —2] the output models have a large spread at these heights. We conclude that 
the inversion results at this height range are not reliable, and that the retrieved atmospheric 
values are reliable only in the fairly narrow height ranges. 

In figure [T] are also shown the retrieved magnetic field strengths with nodes in both the 
field strength (3 nodes) and inclination (4 nodes). The spread of field strengths obtained is 
unreasonably large, especially in the chromosphere. Since the inversion code has no informa- 
tion of Stokes Q and U , some ambiguity between the magnetic field strength and inclination 
cannot be removed. Due to this, for the observations we choose to force the inclination to 
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zero and consider only the longitudinal magnetic flux density which is better constrained. 



2.2. Velocity and inclination discontinuities 

Gradients and/or discontinuities in velocity and magnetic field cause asymmetries in 
Stokes V profiles. The high-/? simulation (P06) demonstrated how velocity gradients from 
shocking acoustic waves cause a pattern of time-varying Stokes V asymmetries in the Ca 
lines. 

To study the effect of discontinuities we use the inversion code in synthesis mode to pro- 
duce a set of Stokes profiles from atmospheres with velocity discontinuities of different sizes 
(1, 3 and 5 km/s) and at different locations (between log(r) = [—5.5,0]. W e then compute 



Stokes V asymmetries for the profiles. The asymmetries are defined as in iMartfnez Fillet 



( 119971 ). The amplitude asymmetry of a Stokes V profile is given by 



dl) -\- CLf. 

where ab and are the unsigned extrema of the blue and red lobes of the Stokes V profile. 
The area asymmetry of a Stokes V profile is given by: 



OA = S—r^ , 

where s is the sign of the blue lobe. 

In figure [2] the results are shown for the velocity discontinuities. For the sign of the 
velocity we use the standard astrophysical notation where positive values indicate red-shifts. 
The sign and amplitude of the Stokes V asymmetries depend both on the location and size 
of the discontinuity. Note that both positive and negative asymmetries are present even 
though all velocities are positive and in all cases the photosphere has a smaller velocity. The 
overall patterns in the 8498 A and 8542 A lines are not identical, or even all that similar. 

To see if the inversion code can handle velocity discontinuities we invert the profiles 
(Stokes I and V of the two Ca lines and the two photospheric Fe lines) using either 3 or 4 
nodes in velocity. Asymmetries computed from the inversion fits are shown in dashed and 
dotted lines in figure [21 In most cases the fits do not reproduce the asymmetries. However, 



-8 -6 -4 -2 2 -a -6 -4-2 2 -fl -6 -4-2 2 
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Fig. 1. — Sensitivity of the lines in the inversions. Upper row: input Stokes / and V profiles 
are in red and fits to the profiles in black lines. Lower row: Temperature (left), longitudinal 
magnetic flux (middle) and velocity (right) as a function of optical depth. Black points show 
results for models with variation in magnetic field strength only, orange points for models 
with variation in both inclination and magnetic field strength. Blue and orange curves show 
the mean of the black and orange points, respectively, and the error bars show the standard 
deviation. 
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Fig. 2. — AsjTiimetries in Stokes V profiles caused by discontinuities in velocity. Upper 
row: Stokes V amplitude asymmetry in 8498 A (left) and 8542 A (right) as a function of 
the location of the discontinuity. Green, blue and red solid lines are discontinuities of 1, 3 
and 5 km/s, respectively. Dotted and dashed lines are asymmetries from inversion fits for 
3 and 4 velocity nodes, respectively. Color coding is the same as for the solid lines. Where 
no inversion value is plotted, no satisfactory fit was found. Lower row: As upper row except 
for Stokes V area asymmetry. 
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the overall qualitative shape of the line profiles is in most cases satisfactory. Increasing the 
number of velocity nodes to 10 does not significantly improve the results. 

Also discontinuities in inclination cause complex asymmetry patterns in the lines. The 
effects of a discontinuity on Stokes V asymmetries are not easily predictable since they 
depend strongly on both the size and location of the discontinuity. Furthermore, the dis- 
continuities can cause the two Ca lines, with relatively close formation heights and similar 
radiative transfer, to exhibit very different asymmetries. Also the area and amplitude asym- 
metries can be of opposite signs. A combination of discontinuities or strong gradients in 
velocity and inclination will be even more difficult to disentangle. 



3. Inversions of quiet Sun observations 

We use the non-LTE inversion code to invert quiet Sun observations of the Ca II 8498 A 
and 8542 A lines and Fe I 8497 A and 8538 A lines taken with the the Spectro-Polarimeter 



for INfrared and Optical Regions (SPINOR. ISocas-Navarro et al.ll2006l ) at the Dunn Solar 



telescope, Sacramento Peak Observatory. The raster scan was obtained in a quiet Sun region 
near disk center at S17.3 W32.1 on May 19, 2005 at 14:14 UT. The upper part of the slit was 
positioned above internetwork and lower part above an enhanced network patch consisting 
of a single polarity. For a detailed description of the observations see P07. 

The inversion is performed in 2 cycles. In the first cycle both Stokes I and V are weighted 
equally and there are 4, 1 and 3 nodes in temperature, magnetic field strength and velocity, 
respectively. In the second cycle more weight (factor 500) is put on the Stokes V profile and 
the number of nodes are increased to 7, 3 and 4 (temperature, magnetic field strength and 
velocity). A node for macroturbulence is also added. The second cycle allows the code to 
fine-tune the solution obtained in cycle 1. In both cycles the magnetic field inclination and 
azimuth are forced to zero. 5 different initial models are used for each inversion and the best 
fit is chosen based on the smallest value. 

A 2-component atmosphere is needed to fit most of the line profiles satisfactorily. The 
filling factor is a free parameter and one of the atmospheric components is a non-magnetic 
external atmosphere which is the same for all pixels. The non-magnetic external atmosphere 
is obtained by inverting an average non-magnetic profile using the same number of nodes 
in velocity and temperature as in the main inversions, but with no nodes in magnetic field 
and a zero magnetic field in the initial models. The filling factor is forced to be unity. The 
non-magnetic profile is an average of intensity profiles in all pixels where the 8498 A Stokes 
V amplitude is below 7 x 10"'^. Because of atomic polarization. Stokes Q and U are given no 
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weight in the inversions. Consequently, we do not have enough information to resolve both 
atmospheric components, hence the non-magnetic atmosphere is set to be constant. 

Since the telluric lines present are weak, they cannot be used for a precise absolute 
wavelength calibration. Instead, the velocities are calibrated so that the mean photospheric 
velocity in "non-magnetic" regions (pixels where the photospheric magnetic field strength is 
less than 300 G) is set to zero. 

3.1. Goodness of fits 

■y^ was used as a proxy for the goodness of the fits produced by the inversions (figure [3]). 
In general, the x^-values are below 0.2. Figure H] illustrates how well the inversions retrieve 
quantities used to describe the line profiles. The inversion code can retrieve the amplitudes 
of both Stokes / (integrated intensity in a 0.75 A wide wavelength band around the line 
core zero wavelength normalized to the band width, /c/A) and V very well, but the Stokes 
V asymmetries are not reproduced. The 8498 A area asymmetry is especially problematic. 
Since the two Ca lines have very different statistics in terms of the area asymmetries (P07) 
and, as shown in the previous section, discontinuities in velocity and/or inclination cause 
complex asymmetry patterns, the failure is not all that surprising. The only quantity that 
shows a correlation with is Stokes V amplitude. This may be explained by it being easier 
to fit only the intensity profiles instead of fitting both the Stokes / and V profiles. Also it is 
possible, if not likely, that the magnetic elements are sub-resolution and the pixel consists of 
a magnetic component and an external non-magnetic component. It is uncertain how well 
the average external non- magnetic atmosphere describes the non-magnetic component. No 
clear correlation is found between and Stokes V asymmetries. 

3.2. Spatial patterns 

Since the inversions have the filling factor as a free parameter, the results need to be 
weighted in respect to how much of the light from the pixel originates from the magnetic 
component and how much from the external non-magnetic component. Since this atmosphere 
is assumed the same for all pixels and only the filling factor varies, for the following discussion, 
only pixels with a filling factor larger than 15 % are considered without further restrictions. 
For pixels with filling factors less than 15 %, only those with 8498 A Stokes V amplitude 
above 0.02 l^. are included. The Stokes V amplitude threshold was chosen to be the amplitude 
above which the standard deviation of in models producing good fits in a given pixel is 
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Fig. 3. — of the inversion fits, right: image of the spatial distribution of logiQ{x^). The 
black strip covers the hairline and the close-by pixels that are affected by it. Left: histogram 
of the values. 
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Fig. 4. — Scatter plots of the observed and fitted 8498 A and 8542 A Stokes / and V 
amplitudes, and Stokes V area and amplitude asymmetries. The over-plotted diamonds and 
asterisks show the mean for 8498 A and 8542 A lines, respectively. In the bottom right 
corners of the figures are cross correlation coefficients for the observed and fitted values. 
Darker gray is 8498 A and lighter gray 8542 A. 
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the same for pixels with filhng factors less than 15% and pixels with filling factors greater 
than 15 %. For pixels with filling factors less than 15 % and 8498 A Stokes V amplitude 
below 0.02, the standard deviation of the inverted atmospheres is significantly larger, and 
they are omitted from the analysis. 

Note that the filling factors are largest in the network and the very quiet internetwork 
whereas the areas surrounding the network have very small filling factors (figure [5]) . This 
illustrates two points: the network is best described with a 2-component scenario. Second, 
since the internetwork has no signal in Stokes V the code can combine the non-magnetic 
profile with a significant contribution from an additional (essentially also non-magnetic) 
atmospheric component resulting in a 2-component non-magnetic atmosphere. In the regions 
by the network, the Stokes V amplitude is in general small, but not negligible. Since the 
Stokes V signal is weighted more by the inversion code in such cases the non-magnetic 
component is described solely by the external non-magnetic atmosphere even though it is 
the dominant component. Since there is little signal in Stokes V in these pixels, consequently 
also the filling factors will be small. 

The photospheric longitudinal magnetic flux density (at log(r)=-l, hereafter referred 
to as the magnetic flux, panel c) in figure [6]) has a fragmented pattern with distinct edges. 
The strongest flux, about 1200 Mx/cm^, is in the centers of flux "bundles". There are no 
isolated points of strong flux. This is due to both the photospheric lines being weak and thus 
not very sensitive to magnetic fields, and to the spatial resolution of the observations. Since 
individual flux tube size is below the spatial resolution, only regions with a concentration 
of several flux tubes appear strong. The strongest flux in the chromosphere (at log(r)=-5, 
panel f) in figure [6]) does not lie directly above the photospheric flux concentrations, and 
not all of the strong flux seen in the photosphere reaches the chromosphere. In fact, the 
strongest flux concentrations are significantly weaker higher up, and the total area covered 
by strong flux is smaller. The magnetic network is more confined in the photosphere. In the 
chromosphere the field has expanded to fill more space causing the more diffuse appearance 
of the network. 

The photospheric temperature (at log(r)=0, panel a) in figure E]) is significantly cooler 
in the strong fiux concentrations. To maintain a pressure balance between the interior and 
exterior of the fiux tube, the gas pressure inside the tube must be smaller than outside leading 
to a reduced density inside the fiux tube. Also the convective energy transport in the tube is 
suppressed by the magnetic field. If the diameter of the fiux tube (or concentration of tubes) 
is large enough the lateral radiative fiux from the hotter surroundings does not reach the tube 
center, and it remains cooler (and appears as a dark feature in intensity images). Hints of a 
granulation pattern are visible in photospheric temperature: the cell centers are hotter than 
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Fig. 5. — Right: filling factors for the inversions. The contours are for 75 G photospheric 
longitudinal magnetic flux. Left: mask based on chromospheric temperature used to define 
the enhanced network, low-flux network and internetwork. 
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the edges. The network is well visible in the chromospheric temperature (at log(r)=-5, panel 
d) in figure [6]), and is clearly hotter than the internetwork. The locations of the hottest pixels 
are not correlated with the strongest fiux. Instead the hot pixels form a fragmented pattern 
that covers most of the region where fiux is seen. The network boundaries are visible as a 
region of intermediate temperature, clearly distinct from both the network and internetwork. 
The intermediate temperature region extends to form a network like pattern in the mostly 
fiux-free region (hence hereafter this region will be referred to as low-fiux network). There 
is no clear indication of the low-fiux network in photospheric temperature or magnetic fiux. 

Since no absolute wavelength calibration was done, defining the zero velocity is not exact 
and velocities close to zero should not be taken as definite measures. In the photospheric 
internetwork patches of negative and positive velocities are present (at log(r)=-l, panel b) 
in figure [H]). The network is dominated by negative velocities. Some positive velocities are 
also seen, though they are located closer to the edges and are not as large in amplitude 
as in the internetwork. There are large scale patterns of positive velocities, interlaced with 
significantly smaller patches of negative velocities, in the chromospheric internetwork (at 
log(r)=-5.5, panel e) in figure [6]). The chromospheric network has mostly positive velocities, 
though above one of the strong fiux concentrations a patch with a large negative velocity is 
seen. The overall velocity pattern in the chromosphere is more fragmented in the network, 
and also the velocity amplitudes vary less than in the internetwork. This may be due to 
the underlying magnetic field "compartmentalizing" the plasma. There is little correlation 
between photospheric and chromospheric velocities. 



3.3. Histograms of atmospheres 

In figure [7| are shown histograms of the photospheric and chromospheric temperature, 
magnetic fiux and velocity. The photospheric temperature histogram is narrower than the 
chromospheric one. It has a single peak at 6400 K, which is also the temperature in the 
photosphere of the external non-magnetic atmosphere. The chromospheric temperature his- 
togram has 2 peaks: the first one is at 5800 K, and corresponds to the low-fiux network. The 
coolest pixels, below 5400 K, are in the internetwork. The second peak, corresponding to the 
hot network pixels, at 6400 K, is narrower but has an extended tail to high temperatures. 
This is higher than in the external non-magnetic atmosphere (6260 K). 

The photospheric velocity histogram is centered at km/s, which is also by definition 
the photospheric velocity in the external non-magnetic atmosphere and in pixels with a 
photospheric magnetic fiux less than 300 G. The chromospheric velocity histogram peaks at 
a positive values: 0.5 km/s. It is significantly wider than the photospheric histogram. 
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Fig. 6. — a) Photospheric temperature (at log((r))=0). b) Photospheric velocity (at 
log((r))=-l). c) Photospheric longitudinal magnetic flux, (at log((r))=-l) d)-f) as a)-c) 
but for the chromosphere (at log((r))=-5, -5.5 and -5). For explanation of choice of pixels, 
see text. Numbers indicate locations of pixels discussed in more detail in section [331 



- 16 - 



The longitudinal magnetic flux in both the photosphere and chromosphere has expo- 
nentially decaying histograms. The flux in the photosphere decays slower than in the chro- 
mosphere. Neither histogram has peaks corresponding to the network flux. 

There are three distinct regions in the chromospheric temperature: the enhanced net- 
work, low-flux network and internetwork. We deflne a mask based on the temperature: the 
internetwork is where the chromospheric temperature is below 5400 K, low-flux network 
where the temperature is between 5400 K and 6200 K and network where the temperature is 
above 6200 K. The mask clearly separates the map into three distinct regions (shown in flg- 
ures[5]and[8]where the three colors mark the three regions). Histograms illustrate similarities 
and differences between the three regions and the line proflles in them. The photospheric 
temperature distributions are nearly identical in all three regions: the network, except for 
the cool strong flux concentrations, is not visible in photospheric temperature. 

The photospheric velocity histograms are centered around zero or close to zero. The 
chromospheric histograms peak at positive values with the internetwork having the largest 
red-shift: it peaks at 1 km/s. The enhanced network regions have the narrowest velocity 
histograms both in the photosphere and in the chromosphere. This is consistent with the 
magnetic fleld suppressing convection. However, the bottoms of the histograms in the en- 
hanced network are roughly as wide as in the internetwork histograms. There is no signiflcant 
difference between the internetwork and low-flux network histograms in the photosphere. In 
the chromosphere, however, the low-flux network histogram is signiflcantly wider. 

As expected, the magnetic flux is by far largest in the network. No flux is seen in the 
internetwork. (Since the photospheric lines are weak and the seeing conditions during the 
observing sequence were not optimal, no internetwork flelds are seen.) The low-flux network 
has some, though very little when compared to the network, magnetic flux. 

The network is clearly brightest of the three regions when viewed in continuum intensity. 
No difference is seen between the internetwork and low-flux network. In the Ca line intensities 
the low-flux network is somewhat brighter than the internetwork, but still clearly darker 
than the network. The Stokes V amplitudes of the Ca lines further illustrate that the 
low-flux network is magnetic: the amplitudes are larger in the low-flux region than in the 
internetwork. In general, the differences between the three regions, especially between the 
internetwork and low-flux network, become more pronounced with height. 

The dependence of the Ca lines' Stokes I amplitudes on the chromospheric temperature 
is shown in flgure M Both lines show a saturation in intensity at high temperatures, and 
a minimum level intensity at low temperatures. The intensities show a strong dependence 
on the temperature only in the intermediate temperature range. Implications on the use 
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Fig. 7. — Left: Histograms of photospheric (dark gray line) and chromospheric (light gray 
line) temperature (left), velocity (middle) and longitudinal magnetic flux (right). 
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Fig. 8. — Histograms of photospheric and chromospheric temperatures, velocities and mag- 
netic flux, continuum intensities and Ca hne amphtudes in the internetwork, low-flux network 
and network. Image shows the mask used to identify the pixels: dark gray is enhanced net- 
work, light gray low-flux network and white internetwork. The black regions are pixels not 
included in the analysis. 
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Fig. 9.— 8498 A and 8542 A Stokes I amplitudes versus chromospheric temperature. The 
network, low-flux network and internetwork are shown in different shades of gray: the dark 
is the internetwork, medium dark the low-flux network and light the enhanced network. 
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of proxies based on spectral lines formed in non-LTE are discussed in section |H Scatter 
plots of chromospheric temperature and longitudinal magnetic flux, both photospheric and 
chromospheric, show some signs of saturation at high fluxes, but the spread is large and the 
number of data points with large magnetic flux is not high enough for the saturation to be 
determined as statistically significant. 

3.4. Individual atmospheres 

In the following, individual atmospheres representative of different regions seen in the 
map are discussed. Locations of the regions are marked in figure El 

Line profiles and atmospheres for pixels belonging to area 1 are shown in figure [TOl Area 
1 is located in a strong photospheric flux-bundle. It is also clearly visible in photospheric 
temperature. The line profiles are all fairly symmetric and similar to one another (except 
for one asymmetric 8498 A Stokes V profile. It corresponds to the atmosphere that deviates 
most from the mean in both the magnetic field and velocity). The average filling factor 
of the pixels is 0.66. The Stokes V profiles of the photospheric iron lines are identical to 
one another indicating that the photospheric magnetic field is fairly homogeneous in area 
1. Photospheric temperature in the atmosphere is 430 K cooler than in the external non- 
magnetic photosphere. The chromosphere is significantly hotter, by 1000 K. The longitudinal 
magnetic flux is of the order of 1 kG in the photosphere and decays by a factor of 5 before 
reaching the chromosphere. This is already seen in the images of magnetic flux where the 
strong photospheric flux concentrations are not visible in the chromosphere. The velocity is 
zero in the photosphere and positive in the chromosphere. 

Area 2 (figure [TTl) is located deep in the internetwork. The Stokes / profiles are all 
nearly identical: deep, narrow and very symmetric. There is no signal above noise in Stokes 
V, and consequently, the filling factor (0.51) only describes how much of the average non- 
magnetic profile is present in the line profiles. The temperature in the photosphere is nearly 
identical to the external non-magnetic atmosphere, but the chromosphere is significantly 
cooler. Because the pixels for the non-magnetic profile were chosen based on magnetic signal 
and not on intensity, pixels in the low-flux network (which are hotter than the internetwork) 
were included. The inversion code returns a flux that is less than 20 Mx/cm^, and the fitted 
Stokes V profiles are well below noise level. The velocity profile differs from the non-magnetic 
both in the photosphere and chromosphere. This is expected since the non-magnetic profile 
is an average of many pixels with different phases, i.e., all but net shifts should be eliminated 
from it. In general, there is very little scatter between the model atmospheres indicating 
that the region is fairly homogeneous, both in terms of thermal and dynamic properties. 
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Fig. 10. — Upper row Stokes / and V profiles for area 1 (for location see figure El). Ob- 
served profiles and inversion fits are shown in black and red lines, respectively. Lower row: 
Temperature (left), longitudinal magnetic flux (middle) and velocity (right) as a function 
of optical depth. The dashed red lines in the temperature and velocity panels are for the 
external non-magnetic atmosphere. The vertical dotted lines show where the inversion code 
results are most reliable. 
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Pixels in area 3 (figure [T2|) are located in the network, but outside of the strong flux 
concentrations. The Stokes / and V profiles are fairly symmetric, but only the 8542 A 
line has nearly identical Stokes V profiles in all pixels. The other 3 lines show two sets of 
Stokes V amplitudes: one set significantly larger than the other. This is especially clear 
in the photospheric iron lines. The 8498 A line has a set of Stokes V profiles with strong 
signal in the line wings and a set with little signal in the wings and also a slightly smaller 
amplitude. The temperature in the photosphere is roughly the same as in the external non- 
magnetic atmosphere: the coolest photosphere in area 3 is only 50 K cooler. This is a very 
small difference when compared to area 1 where the Stokes V profiles had similar shapes 
(albeit somewhat larger amplitudes) as here. In contrast to area 1 where the chromosphere is 
significantly hotter than the external non-magnetic atmosphere, in area 3 the chromosphere 
is cooler, by roughly 150 K. Above log(r) = —6 the temperature decreases even further. This 
is likely a consequence of spline interpolation when the code is adjusting the temperature 
in the lower parts of the atmosphere. Since the lines are not sensitive to conditions above 
logir) = —6 this peculiar temperature profile has no effect on the fits. As expected based on 
the two sets of Stokes V amplitudes, there are two sets of magnetic fluxes in the photosphere. 
The spread of fluxes in the chromosphere is significantly smaller. Both the photospheric and 
chromospheric velocities are lower than in area 1. The velocities are slightly negative in the 
photosphere and positive in the chromosphere. The models are in fairly good agreement with 
one another. Since the code cannot retrieve sizes of velocity gradients, this does not rule 
out that the differences in the 8498 A and 8542 A profile shapes are not caused by large 
gradients or discontinuities in velocity (or inclination). 

Area 4 (figure [T3!) is also located in the network, but even farther from the strong 
photospheric magnetic flux concentrations than area 3. The chromospheric Stokes I profiles 
show self- reversals which are not seen in area 3. The self-reversals are stronger on the red 
side in the 8542 A line and in most cases also in the 8498 A line. The chromospheric 
Stokes V profiles all show varying degrees of asymmetry where the red sides tend to have 
multiple components. The inversions reproduce the general shapes of the profiles, but not the 
amplitudes of the self-reversals and the emission features in Stokes V. Since the external non- 
magnetic atmosphere is based on an average profile of both the map and time series, anything 
leading to asymmetries in Stokes / has to be incorporated into the magnetic component. The 
photospheric temperature is not very different from the external non-magnetic atmosphere 
but the chromosphere is clearly hotter. 
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Fig. 12. — As figure [TOl but for 
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4. Discussion 

Two points need to be kept in mind when interpreting the inversion results. Firstly, the 
sensitivity of the inversions is limited to specific heights. Since the Fe lines and Ca line wings 
are formed in the photosphere, and only the Ca line cores are chromospheric, the sensitivity 
range is larger in the photosphere than in the chromosphere. Because the inversion code 
does not include treatment of atomic polarization (and consequently Stokes Q and U cannot 
be used in the weak field regime), it is not possible to resolve the ambiguity between the 
magnetic field strength and inclination, instead only the longitudinal magnetic flux density 
can be retrieved. 

Secondly, discontinuities in velocity and magnetic field inclination result in complicated 
Stokes V asymmetry patterns. The responses to discontinuities in the two Ca lines are not 
always similar. Also the area and amplitude asymmetries can respond in different ways. The 
inversion code cannot reproduce the discontinuities. Such discontinuities may explain the 
observed Stokes V area and amplitude asymmetry histograms (P07): the 8542 A area asym- 
metry histogram peaks at a positive value whereas all other histograms (8542 A amplitude 
asymmetry, 8498 A area and amplitude asymmetry) peak at a negative value. Based on the 
high-/5 simulations it seems unlikely that the difference is caused by shocking acoustic waves 
alone. The current inversions are not a good diagnostic for discontinuities or strong local- 
ized gradients, especially since both velocity and inclination discontinuities can be present. 
Since we do not include linear polarization in the inversions, and since the asymmetries are 
not reproduced by the inversions, there is information in the spectra which is not used by 
the current inversion code. A more suitable inversion code, especifically designed for the 
inhomogeneities of the chromosphere, is needed to solve the problem. 

Despite the limitations, these non-LTE inversions remain a useful tool: observed Stokes 
/ and V amplitudes are reproduced successfully even if asymmetries in Stokes V are not. 



4.1. Magnetic field 

The magnetic contrast of the network and surroundings is significantly reduced in the 
chromosphere, and there are clear indications of the network expanding with height. These 
indications include the spatial distribution of Stokes V amplitudes in the Ca and Fe lines, the 
appearance of the magnetograms (figure 2. in P07) and spatial distributions of temperature 
and magnetic flux in the photosphere and chromosphere. The dif fusive appearan c e and 



spreading of the network in the chromosphere was already noted by iHarvey &: Hall! (jl97ll ) 



who attributed it to the magnetic field becoming more horizontal with height. This is in 
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agreement with our data, and especially with the clearly more diffuse and larger appearance 
of the network in the 8542 A magnetogram when compared to the 8498 A magnetogram. 
However, since we do not have information of the field inclination, we cannot conclude this 
for sure. 

Not all of the strong magnetic flux concentrations in the photosphere penetrate into 
the chromosphere. One possible cause is the expansion of the field with height. This would 
reduce the longitudinal flux and increase the amount of horizontal fleld on the expense of 
vertical fleld. Another e x planation is that s o me of the flux returns bac k to th e photosphere. 
Schrijver fc Title! (120031 ) , IJendersie fc Peterl (120061 ) and lAiouaz fc RastI (120061 ) have explored 
the connection of the network magnetic fleld to the corona, and shown that the expansion of 
the network is affected by magnetic flelds of opposite polarity in the network and internet- 
work. Unfortunately, our observations of the photospheric magnetic fleld are greatly limited 
by the weakness of the iron lines. There is no indication of opposite polarity magnetic fleld 
in our data, but as pointed out by the above-mentioned authors, sub-resolution magnetic 
flelds in the internetwork play a crucial role in determining what fraction of the network flux 
survives to the corona and what fraction returns to the photosphere forming loops at various 
heights. The low-flux network is a clear indication of there being a substantial amount of 
magnetic flux not captured by the photospheric iron lines. 

Based on the inspection of the individual atmospheres the photosphere appears more 
homogeneous, especially in terms of magnetic flux. There is more small-scale spatial variation 
in the chromosphere. The photosphere is unlikely to have large gradients or discontinuities, 
and the photospheric radiation is in LTE. Consequently, it is easier to flnd a model to flt 
the observations in the photosphere and the spread of atmospheric models will be smaller 
there. However, since the small-scale chromospheric spatial patterns are coherent over more 
than just one pixel, and because we see spatial variations on small scales in the line proflles 
as well, it is likely that the inhomogeneous appearance of the chromosphere is real. 



4.2. Dynamics 



Velocities in the chromosphere, especially in most of the network, are predominantly 
red-shifted in relation to the photosphere. The chromospheric network velocity patterns are 
fragmented and large velocities are more often seen outside the network. The latter is also 
true in the ph otosphere wh e re the difference is even more pronounced. A similar flnding 
is presented in iRezaei et al.l (120071 ) where the fraction of large velocities (determined from 
photospheric Stokes V proflles) is larger in the internetwork than in the network. In the 
current data, there is no difference between the low-flux network and internetwork in terms 
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of photospheric velocities, but in the chromosphere the distribution of velocities is slightly 
wider in the low-flux regions. The strongest magnetic fields suppress oscillations, although 
large asymmetries are also seen in network line profiles. In general, the largest asymmetries 
are located at network boundaries. Also increasing Stokes V asymmetries are correlated 
with decreasing magnetic flux. The asymmetries are in most cases at least partly due to 
emission features in the line profiles, i.e. they are also likely to be related to the temperature 
structure, not just velocity and magnetic field gradients. The difference between the velocity 
histograms in the low-fiux network and internetwork indicates that magnetic fields that do 
not play a visibly important role in the photosphere, become important higher up. 



Recent work by e.g., IJefferies et all (120061 ) and lVecchio et al.l ( 120071 ). has tied together 
the magnetic topology and the observed oscillatory properties in the chron aosphere. Th e 



network edges being more dynamic is consistent with the picture presented by IJefferies et al. 



(120061). In it i ncline d magnetic fields channel low-frequency waves into the chromosphere. 
Jefferies et al.l ( l2006l ) also point out that the low-frequency waves are intermittent owing to 
the continual changes caused by the motion of the convective cells and by the "magnetic 
carpet" (i.e., small magnetic dipoles created in the internetwork and transported to the 
network where they interact with existing flux). Reconnection events would also probably 
lead to asymmetries in line profiles. In our data the network clearly expands with height, 
and the difference between the 8498 A and 8542 A magnetograms is striking, indicating that 
conditions described above (i.e., inclined fields and more dynamic network boundaries) are 
likely to exist. Furthermore, since the photospheric driver has most of its power at periods 
of 5 min, the driver of chromospheric oscillations in regions with inclined fields is stronger 
than in the internetwork where the 5 minute period waves are filtered out. This was pointed 
out by iDe Pontieu et al.l (120071 ) in connection to the dynamic behavior of fibrils. All the 
above is consistent with the network, and especially the boundaries, being more dynamic 
and exhibiting more self-reversals and profile asymmetries than the internetwork. 



4.3. Temperature proxies and heating of the magnetic chromosphere 



Figure [9] shows that the intensity of the Ca II IR triplet lines saturates at high chro- 
mospheric temperatures. A si milar saturat i on eff ect is seen if the magnetic flux is plotted 
instead of the temperature. In lRezaei et al.l (120071 ) a similar plot, but for the Ca H intensity 
and photospheric flux, is interpreted as a saturation limit for the magnetic heating (i.e., the 
heating as a function of field strength or filling factor levels off at high values). The inver- 
sions show that the leveling off in the Ca II IR triplet lines is, at least partly, a radiative 
transfer effect, and does not necessarily imply that the magnetic heating has an upper limit 
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that depends on the field strength or fiUing factor. The data do not conclusively show this 
not to be the case either. (In fact, there are some signs of a saturation effect in scatter 
plots of magnetic flux and chro mospheric temperatu re.) Also a lower limit on chromospheric 



emission, similar to that seen in lRezaei et al.l (120071 ). is seen in our data, but no lower limit 



on the chromospheric temperature is found. The Ca II IR triplet line intensities are linearly 
dependent on temperature only in the range between roughly 5800 and 6500 K. Interpret- 
ing proxies based on lines formed outside of LTE is not straight-forward, and the effects of 
radiative transfer need to be considered in detail when inferring physical properties, such as 
velocity and temperature, from proxies. 

The network patch appears on large scales fairly uniform in temperature, and the hottest 
pixels do not always coincide with strongest photospheric or chromospheric flux. Nor do they 
coincide with the largest fllling factors. Based on this the relationship between heating and 
magnetic flux does not appear to be linear. However, since the low-flux network is cooler, 
the amount of flux is clearly a factor in the heating. We see signs of a satura tion effect in 



the he ating (flux as a function of chromospheric temperature) as proposed by iRezaei et al. 



(120071 ) and references therein. The number of data points is not, however, large enough for 
making deflnite conclusions. Further observations and inversions, preferably simultaneously 
with strong photospheric lines, are needed. 



5. Final remarks 

Because of non-LTE it is important to take into account radiative transfer effects when 
interpreting proxies based on line proflles. This is why non-LTE inversions are an invaluable 
tool: with them it is possible to infer macroscopic physical quantities without resorting 
to proxies. The inversions presented in this work further highlight the physical differences 
between the dynamics and structure of the chromospheric network and internetwork. The 
chromosphere is divided into three regions characterized best by chromospheric temperature, 
but also visible in other chromospheric parameters. The fundamental differences between 
the three regions are most likely related to magnetic flelds. 

There are limitations to what can be done with current non-LTE inversion codes, and 
forward modeling remains a crucial part in understanding the dynamics of the magnetic 
chromosphere. We plan to make simultaneous observations of the Ca II IR triplet lines and 
the He I 10830 A lines to further validate the inversions, to obtain a more 3-dimensional 
view of the magnetic chromosphere, and to study the connection between thermal properties 
of the chromosphere (magnetic heating) and magnetic fleld topology. 
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